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Abstract. - Using the recently developed so-called White Bear version of Rosenfeld's Funda- 
mental Measure Theory we calculate the depletion potentials between a hard-sphere colloidal 
particle in a solvent of small hard spheres and simple models of geometrically structured sub- 
strates: a right-angled wedge or edge. In the wedge geometry, there is a strong attraction 
beyond the corresponding one near a planar wall that significantly influences the structure of 
colloidal suspensions in wedges. In accordance with an experimental study, for the edge ge- 
ometry we find a free energy barrier of the order of several ksT which repels a big colloidal 
particle from the edge. 



When a colloidal suspension contains particles of different size, depletion forces arise be- 
tween big particles due to excluded volume effects. If a colloidal suspension is prepared such 
that electrostatic interactions are screened and only hard-core repulsive interactions are left, 
these forces are attractive at small distances pQ and can exhibit a potential barrier and an 
oscillatory tail for large separations PUSSHS] • The same forces also arise if a big colloidal 
particle b (diameter er&) suspended in a solvent of small particles s (diameter a s ) approaches 
a wall. The depletion forces for hard-sphere models of bulk colloidal suspensions, as well as 
colloidal suspensions close to structureless walls have been the subject of many theoretical 
studies. Technical developments during the last decade made it possible to prepare well- 
defined, monodisperse colloidal suspensions and with the help of video-microscopy [Sj and 
total internal reflection microscopy [7|03EI1 such depletion potentials can be now measured 
directly and compared quantitatively with theoretical predictions. 

On the other hand, the behavior of colloidal suspensions close to structured walls and 
resulting depletion potentials have not yet been investigated in depth. Such situations arise 
naturally for example close to a corner in a container for colloidal suspensions. Geometrically 
structured substrates with a dedicated design have been used to control the growth of colloidal 
crystals ^U] . In this context Dinsmore et al. studied the motion of colloidal " hard-spheres " 
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near the edge of a terrace by means of video-microscopy. In the direction perpendicular to 
the edge they found a free-energy barrier approximately two times the mean thermal energy. 
Dinsmore ct al. suggested that localized cntropic force-fields created by various structures of 
the substrate can be used to control the movement of colloidal particles and, consequently, 
could be used to create highly ordered arrays of colloidal particles. 

Recently, Kinoshita et al. ^JEj have studied depletion potentials of hard spheres close 
to the substrate exhibiting various geometrical features by means of first-order Ornstein- 
Zcrnike (FOOZ) integral equations. Using correlation functions of the bulk hard-sphere fluid 
as input, they calculated the potential of mean force between a big particle and an edge or 
hemispherical hole in a planar substrate (mimicking key and lock steric interactions between 
macromolcculcs). However, FOOZ theory assumes that the correlation functions close to the 
surface are the same as in the bulk, which gives rise to discrepancies in predicting packing 
effects as compared to simulation data. For geometrically structured surfaces, where these 
packing effects may be very pronounced like, e.g. in a wedge, the FOOZ approach is expected 
to be even less reliable. 

Depletion forces in the bulk or close to planar walls have been studied also by means of 
molecular dynamics [2] and Monte Carlo (MC) simulations While being computationally 
quite demanding, these studies serve as a useful test for various theories. However, the study 
of the fluid structure near geometrically structured substrates is computationally even more 
costly than for planar walls so that there are only few such studies |14lll5llTO| . 

A versatile approach to study depletion potentials close to arbitrarily shaped substrates 
has been proposed in Refs. 0] and [5]. According to this theory the depletion potential is 
evaluated using the potential distribution theorem [HI51 ll7llT%] on the basis of the density 
profile of the small spheres close to the substrate without the big sphere. This method is valid 
for all geometries of the substrate, but explicit calculations so far have been carried out only 
for structureless planar walls @] or spherically curved substrates 

Here we apply the aforementioned scheme to calculate the depletion potential of a single 
big hard sphere close to a right-angled wedge or edge. Within this approach it is important to 
determine accurately the density profile p s (r) of the small spheres close to the edge or wedge 
where it varies most rapidly. 

The calculation of the density profile of a hard-sphere fluid close to planar or spherical 
structureless walls constitutes a standard problem in the statistical mechanics of fluids which 
lends itself to be solved by density functional theory (DFT). It has emerged that Rosenfeld's 
Fundamental Measure Theory (FMT) [20] is one of the most accurate DFTs if compared with 
simulations for either one-component, binary [22 or polydisperse [22] hard-sphere fluids close 
to a hard-wall. However, for complex geometries like wedges or edges there are only a few 
comparisons of DFT and simulations available [221E!! ■ Thus before calculating the depletion 
potentials we test the accuracy of the DFT against MC simulations. 

The key approximation in every density functional theory resides in an expression for the 
excess (over the ideal gas) intrinsic free energy J-~ex S \pi(r)] of an inhomogeneous fluid as a 
functional of local densities pi (r) . Within the framework of FMT [20] one has 



where (3 = (ksT) 1 , ks is the Boltzmann constant, and T is the temperature. The excess 
free energy density $ is a function of weighted densities 
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Fig. 1 - Schematic representation of a right-angled wedge z w (x) = The dotted line indicates a 
possible cut along which the density profiles are presented. 

Fig. 2 - Density profile p(x = 0, 2*) from MC (circles), RF version (dashed line) and WB version of 
FMT (solid line) for a bulk density p s = 0.7, with z* — — -j=. The inset illustrates this density 

distribution by showing the corresponding typical particle configuration. Here and below lengths and 
number densities are given in units of a s and a^ 3 , respectively. 



with six different geometrical weight functions w\ per component i (four scalar and two 
vector- like) . In the present problem all quantities are translationally invariant along the 
wedge or edge which we choose as the y-direction. In order to calculate weighted densities 
n a (x, z) in an efficient way we have applied the two-dimensional Fast Fourier Transform. 

From a variety of expressions for 4> here we select two of them. The first one is the original 
Rosenfeld functional (RF) HJ 

$ ( )({»„}) = -no ln(l-« 3 )+ + 24?r(1 _ n3)2 , (3) 

based on the Percus-Yevick compressibility equation of state. As pointed out in Rcf. [25] 
this approximation does not render a freezing transition. Our second choice is based on the 
Boublik-Mansoori-Carnahan-Starling-Leland equation of state as proposed by Tara- 

zona |2H| and Roth et al. (see also Ref. j3U] where the same functional has been proposed, 
too). Moreover this choice implements a modification proposed in Ref. [23] (the so-called 
antisymmetrized form, Q = 3) that allows for the occurrence of a freezing transition: 

ih(WB)fr n w n \ 1 n i n 2 - ni • n 2 3 2 3 n 3 + (1 - ?i 3 ) 2 ln(l - n 3 ) 

® (W >({n a }) = -n ln(l-n 3 ) + +nt(l-£y — — w , 4 

1 — n 3 6birn 3 (l — n 3 ) z 

where ^(r) = |n 2 (r)|/n 2 (r). Following the nomenclature in Ref. this DFT version is 
denoted WB ("White Bear"). The density profiles are obtained by solving the Euler-Lagrangc 
equations corresponding to the minimum of the grand potential functional. 

The Monte Carlo simulations have been performed in a grand canonical ensemble according 
to a procedure described in detail in Ref. ^Hl- However, we emphasize that unlike in Ref. 
the MC-generated density profiles presented below have not been smoothed but are instead the 
"raw" data generated in MC simulations with no further post-simulation processing applied. 
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Fig. 3 - The density profile p(x = 0.36, z*) from MC simulations (circles), RF version (dashed line), 
and WB version (solid lines) evaluated for the bulk density p 3 = 0.7. The inset shows the contact 
value of the density profile resulting from the White Bear version as function of x. 

Fig. 4 - The density profile p(x = 0.5, z*) from MC simulations (circles), RF version (dashed line), 
and W B version (solid line) evaluated for the bulk density p s — 0.7. 



The comparison of the DFT and MC density profiles has been carried out only for the 
wedge geometry, because in this geometry packing effects are much more pronounced than 
for the edge and thus presents the more severe challenge for DFT. The dependence of the 
density profile p s = p s (x,z) on two spatial variables is illustrated by vertical cuts along the 
lines x = const (see fig. QJ. In the following all length scales are expressed in units of a s , the 
diameter of the small spheres, and the number densities in units of <r~ . 

Figures T2I4I show such MC density profiles (circles) and the corresponding ones from DFT 
(dashed and solid lines). The comparison with the MC data shows that the WB version 
represents a noticeable improvement over the RF version in particular for z* = — -1= < 
0.5. Along the line x — (see fig. EJ) the contact value of the profile resulting from the WB 
version is 15.54 while from RF 17.29, i.e., a difference of ca. 10%. Also for x = 0.36 (see fig. 13 
the WB version agrees better with the MC data than RF again in particular for z* < 0.5. 
The contact value p^ B \x = 0.36, z* = 0) = 2.55 differs from p { s RF \x = 0.36, z* = 0) = 1.99 
by ca. 20%. 

The large variations (between 15.5 and 2.5) of the contact densities for < x < 0.5 (see the 
inset in fig. El illustrate the strong packing effects in this system. In fig.EJfor x = 0.5 a small 
plateau develops around z* ~ 0.5. The RF version exaggerates this feature considerably and 
leads to a contact value p s (x = 0.5, .z* = 0) = 1.81 whereas the WB version reproduces 
the plateau more accurately leading to a contact value p^ B '{x = 0.5, z* = 0) = 2.38. 

We refrain from comparing directly contact values from DFT with those from MC because 
in the MC simulations p (x, z) has to be calculated as a histogram on a two-dimensional grid. 
Even though the mesh width Sx = Sz = 0.02 of this grid is small, an extrapolation scheme is 
required to obtain an estimate for the contact value. This extrapolation scheme is particularly 
delicate near the core of the wedge where p (x, z) varies rapidly. However, inspection of figs-EJ- 
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Fig. 5 - The depletion potential of the big colloidal particle close to the edge along the line of 
closest contact (see the inset in fig.©). The packing fraction of the small spheres is 77 = 0.366. The 
dashed, solid, and dashed-dotted lines correspond to the size ratios q = a s /at = 0.333, 0.2, and 0.1, 
respectively. 

Fig. 6 - The depletion potential along the line of closest contact between the big colloidal particle 
and the 90° right-angled edge (inset) evaluated for Oh = 5.35cr s and r\ = 0.3 



0] leads to the expectation that the WB version of DFT reproduces the actual contact values 
rather accurately. We have carried out calculations also for lower bulk densities finding a 
similar (or even better) agreement between DFT and simulations. Thus we conclude that 
FMT predicts accurately the structure of the hard-sphere fluid in a hard wedge; its WB 
version improves the agreement with simulations near the core of the wedge, where packing 
effects are particularly pronounced. 

Reassured by the quality of the density profiles resulting from DFT, we have calculated 
the depletion potentials W near a hard wedge and edge based on the WB version of DFT. 
We use the fact that in the limit of a dilute solution of big particles fiW can be calculated 
from the density distribution of the small spheres unperturbed by the presence of the big one. 
It has been shown 0J[S] that in the dilute limit of the big spheres the depletion potential of 
species & at a point r can be expressed in terms of the difference between the one-particle 
direct correlation function cj; in the bulk (r — > 00) and at r: 

f3W(r)= lim (^(r-.ooj-c^r)) (5) 

with 4 1} (r) = -6F ex /5p b (r) within DFT Ell- 
in fig. [3] we show (3W for size ratios q^ 1 = <J\,ja s = 3, 5, and 10 at a hard edge for a bulk 
packing fraction of the small spheres rj = ^p s cr^ — 0.366. The depletion potential is evaluated 
along the line of contact between the big sphere and the edge (see the inset in fig. 0); L 
parametrizes a path on this line. For L — > ±00 the depletion potential between a planar hard 
wall and a big sphere is recovered. Upon approaching the edge, oscillations in the contact 
value of PW occur. At the corner there is a pronounced maximum in the (still attractive) 
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Fig. 7 - The depletion potential of the big colloidal particle close to the wedge taken along the line 
of closest contact. The packing fraction of the small spheres is r\ — 0.366. The dashed, solid and 
dotted-dashed lines correspond to the size ratios q = 0.333, 0.5 and 0.1, respectively. 

depletion potential. The value of (3W at the maximum, i.e., at L = equals —1.5, —1.6 and 
—2.0 for q = 0.333, 0.2, and 0.1, respectively. This maximum represents an effective repulsive 
barrier repelling a big colloidal particle approaching the edge from the side and practically 
preventing it from passing around the corner. The difference (3AW between the depletion 
potential at the planar hard wall and at the corner is a quantitative measure of the barrier 
height. The barrier height increases upon decreasing q and for the system shown in fig [5] it 
equals 2.1, 3.8, and 8.0, respectively. 

In fig. El we present f3W along the line of contact for a system mimicking the experimental 
system studied by Dinsmore et al. Theoretically we find (3 AW ~ 3.6 which is larger than 
the corresponding result obtained from the Asakura-Oosawa approximation ((3 AW ~ 3, 
which treats the solvent as an ideal gas. The experimental value quoted for the barrier height 
(— 2) JJI is smaller than both theoretical results. This is expected because for the theoretical 
value the big particle is assumed to move along the line of contact with the wall whereas 
experimentally on average some of the big particles monitored experimentally happen to be 
located slightly above the surface. (3 AW reaches a value of 2 at a normal distance of about 
0.2er s . The depletion potentials shown in figs. El and El exhibit an oscillatory behavior due 
to correlation effects missed by the Asakura-Oosawa approximation which holds only in the 
limit r) — > 0. These correlations lead to a minimum of (3W at L ~ 3.5 next to the main peak 
(see fig. El , which is the global minimum and is also visible in the experimental depletion 
potential [TTj . 

In fig. we show (3W along the line of contact for a single big hard sphere in a right-angled 
hard wedge for r\ = 0.366 and q = 0.333, 0.2, and 0.1. As in fig. El we focus on the depletion 
potential along the line of contact between the big sphere and the wedge surface (see the inset 
in fig. UJ. Again in the limit L — > ±oo the contact value of f3W between a planar hard wall 
and a big hard sphere is recovered. Near the core of the wedge oscillations occur which are 
more pronounced than for the corresponding edge cases. Attraction is strongest directly in 
the corner where (3W = —7.5, —10.2, and —20.3 for q = 0.333, 0.2, and 0.1, respectively. This 
is roughly twice the corresponding contact values at a planar hard wall (-3.6, -5.4, and -10.0). 
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This potential is expected to promote highly ordered structures of big colloids in the core of 
the wedge. The pronounced oscillations of the depletion potential upon approaching the core 
of the wedge along the walls can lead to additional ordering. 

We have studied depletion potentials near geometrically structured substrates. It has 
turned out that the recently developed "White Bear" version [2H| of Fundamental Measure 
Theory in its antisymmetrized form |25) is able to reproduce accurately the structure of the 
hard-sphere fluid in a hard wedge as obtained from Monte Carlo simulations, wall. The free 
energy barrier repelling a big colloidal particle from a substrate edge can be of the order 
of several kgT in accordance with experiments In conclusion, the results presented in 

this letter demonstrate that geometrically structured substrates create strong and oscillatory 
entropic force fields that may significantly affect the structure and hence the properties of 
colloidal suspensions near their container walls. 
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